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Abstract 

We present a quantum mechanical study of the diatomic hydro- 
gen H( 2 S) and oxygen 0( 3 P) collision and energy transfer for its four 
molecular symmetry (X 2 U 2 £~, 4 IT, which is important for the 

investigation of many processes of astrophysical and chemical inter- 
ests including the one on molecular cooling, trapping or Bose-Einstein 
condensation. We compute the rovibrational spectra for the (X 2 U) 
state and the resulting bound states are in an excellent agreement 
with experimental data for the energy range lower than the dissoci- 
ation threshold. We calculate the phase shifts of its partial- waves, 
the total cross section, and the differential cross section. They are all 
well-structured because of the shape of its potential curve. We do the 
similar studies for the other three dissociative states ( 2 £ - , 4 II, £"""). 
Finally, we have also decided the thermal rate coefficients for the hy- 
drogen and oxygen collision for its four individual states. 
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1 Introduction 



The atom and atom collisions are the simplest nevertheless among the most 
important systems for quantum molecular scattering investigation. On one 
hand, they serve as a prototype for the understanding of quantum scattering 
process because they are the only real molecular systems at present time 
that a complete picture of the different aspects of the quantum scattering 
can been seen from a computational simulation, and that a strict test with 
the experimental observations can be made |TJI2]- On the other hand, good 
understanding and accurate calculation of related quantities for the diatomic 
scattering are essential to the exploration of many other important fields. 
This includes the study of laser cooling, magnetic trapping and Bose-Einstein 
condensation of molecules [HUE]. To these ends, we have chosen the system 
of oxygen-hydrogen collision for our present study, 

0* + H^O + H*. (1) 

On the practical or application side, this reaction constitutes the basic or 
primary step for many complicated processes for combustion, chemical reac- 
tions, and atmospheric or astrophysical sources [5-30] . For example, it is one 
of the main reactions occurring in the nonthermal escape of the hydrogen 
from Venus and Mars, where O* is produced due to the dissociative recombi- 
nation of O2 in the Martian atmosphere [HI EH EH] ■ Another example is 
that the collisional hydrogen and oxygen has been investigated to illustrate 
the use of infrared laser pulses to manipulate the collision pairs [20] . 

In this paper, our object is to present a relatively detailed and systematic 
study of the collision for the system (1). In Section 2, we first show the 
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potential energy curves of the first four states of the collision (1) and decide 
the bound states for the ground state due to the existence of an attractive 
potential well. We then present the fundamental theory for the investigation 
of their quantum scattering process. In Section 3, we demonstrate the vari- 
ous aspects of the scattering processes. We first show the behavior of most 
fundamental quantity - phase shift - for the quantum scattering and study its 
variation in terms of the partial waves or collision energy. We then calculate 
the total cross sections and the differential cross sections and study their 
relations with the change of energies or scattering angles. We also compute 
the thermal rate coefficients which is the average of the reactive flux over 
the Maxwell-Boltzmann distribution. The final Section contains a discussion 
and summary. 

2 Potential Energy Curves and Scattering The- 
ory 

2.1 Potential energy curves 

We first present the potential energy curves for the collisional 0( 3 P) + H( 2 S) 
system. It is well-known that the dynamic features, especially the resonance, 
is decided by or even sensible to the details of the shape for the curves |3*T] . 
In addition, the previous calculations [T5J [TBI Mj\ for the excited states 
( 2 S~, 4 I1, and 4 S~) use the potential data which are approaching the asymp- 
totic region in the positive direction which is in violation of the physical law 
of an attractive dispersion force at a large distance El 13 El H31 There- 
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fore, we decide to use the one from the calculation of Yarkony without an 
inclusion of any fine-structure corrections such as the spin-orbit interaction 
^3123]. We do the analytical fitting and require that the asymptotical po- 
tentials satisfying the negative or inverse 6th-power of internuclear distance. 
Figure 1 displays the first four states (X 2 U 2 S - , 4 II, 4 £~) of the system (1) 
with the energy zero chosen to be at the asymptotes. The curve for the X 2 U 
state is the typical one of a ground state for a diatomic molecule, featuring 
an attractive potential well at a stable configuration. The other three are 
the usual repulsive potentials seen for the excited states. 

The molecular spectra of the OH system is still of fundamental inter- 
est EH] • In addition to the continuing investigation of its whole range 
of the spectra, the OH maser is considered to be a means for detecting the 
interaction between the supernova remnant and molecular cloud which is im- 
portant for understanding the formation processes of new stars |33J E3 HIE] • 
For the ground state (X 2 U) potential energy curve, we have determined the 
bound states associated with the attractive potential energy well for the en- 
ergy below the dissociation threshold. We do the calculation by utilization 
of the LEVEL 7.5 program [37]. The Hamiltonian includes the standard 
contributions for a diatomic molecule. The resulting rovibrational energy 
levels constitute a banded structure with each vibrational quantum number 
n associated with a bunch of rotational sublevels {j}. The table 1 shows the 
computed results for vibrational quantum number up to 16 and when the j 
equal to or 1. When compared to the experimental observation [USES], 
they have an excellent agreement except for the states near the dissociation 
threshold where we expect a big derivation from the behaviors for the ideal 
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models. 



2.2 Scattering Theory 

We employ the usual approach of partial wave expansion to study the atom- 
atom collisions where both the incident and the scattered waves are analyzed 
in terms of the components of spherical waves |2] . Since only the phase 
shifts 5i of the outgoing partial waves are changed, they are the central quan- 
tity that need to be calculated. They contain all the dynamic information for 
the scattering process. The phase shifts are determined by the requirement 
that the inside and outside scattering waves and their derivatives are contin- 
uous at the boundary point R where the interaction potential V(r) between 
two atoms vanishes. They are given by the formula 



and ji or ni are the Ricatti-Bessel or Ricatti-Neumann functions. The func- 
tion u(R) is the radial scattering wave obtained by numerically solving the 
radial Schrodinger equation from R = for the given central potential. Once 
the phase shifts are known, the differential or the total cross sections, for 
instance, can be easily calculated as follows 



tan Si 



'i = 



u{R)k] l '{kR)-u'(R)] l (kR) 
u{R)kni{kR) - u'{R)ni{kR) ' 



(2) 



where 




dQ 



(3) 



where 



1 



oo 




1=0 
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and 



47T °° 

fc z=o 



(5) 



3 Calculation Results 

3.1 Partial wave analysis 

We first do the partial wave analysis of the ground state (X 2 U) for the O + H 
collisional system. Figure 2 illustrates the variation of the calculated phase 
shifts for the s and p waves in terms of the total collision energy. We observe 
the expected phase change when the angular momentum quantum number 
changes by one unit. For the energies computed, there are many partial 
phases contributing to the total cross section for the scattering process. In 
addition, they have the trend for approaching the zero when seen from the 
lower energy side. 

For the excited states, the partial wave analysis can be carried out in the 
same way. Figure 3 is related to the change of the s and p waves in terms 
of the collision energy for the 4 S~ state. We notice the similar behaviors to 
those for the X 2 Tl state. 

3.2 Total and differential cross sections 

Figure 4 depicts the total cross sections for the four electronic states of the 
hydrogen ( 2 S) and oxygen ( 3 P) collisional system with different symmetries. 
The energy range is from 0.27 eV to 15 eV. As expected, they are a decreas- 
ing function of the total energy except at the points where resonance occurs 
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for the ground state. It is due to the existence of a potential well for this 
state (X 2 I1). In addition, the magnitude of the total cross section for the 
X 2 U state is significant larger than those for the other three repulsive states. 
The latter states show the similar behaviors. 

We also compute the differential cross sections for the collision (1) with 
different orientations, which is a key quantity, for instance, for the study 
of escaping issues of hydrogen and deuterium in the planetary exosphere. 
Figures 5 displays the calculated results for the X 2 U state for the energy 
equal to 0.2 eV or 3.0 eV, respectively. We see that for the low energy there 
is a very large probability for forward scattering (9 = 0) and also a noticeable 
amount for the head-on collision (9 = n). When energy increases, however, 
the chance for the forward scattering persists but the one for the head-on 
collision diminishes. The other states have the similar situations but they 
have no probability for head-on collision for any values of energies. This is 
shown in Figure 6. 

3.3 Reaction rate coefficients 

The thermal reaction rate coefficients at temperature T are defined as follows 



which are the Maxwell-Boltzman average over the reactive flux. The reactive 
flux is the relative velocity v for the collision times the total cross section a 
calculated before. The E is the relative translation energy for the diatomic 



stant k B times the absolute temperature T. The calculated temperature 
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ma mum 





system given by E = ^jj,v 



2 , and the j3 is the inverse of the Boltzman con 



dependence of the reaction rate coefficient for the 4 IT state is demonstrated 
in Figure 7. The temperature range chosen is of planetary interest. We see 
that the calculated result shows an increasing function of the temperature. 
The other states have the similar behavior. 

4 Discussions and Summary 

We have decided the bound states associated with the existence of an attrac- 
tive potential well for the ground state (X 2 n). The computed rovibrational 
energy levels are in an excellent accord with the experimental data for the 
lower energy range. However, when the energy increases up to the direction 
of the dissociation threshold, the difference is accumulating and becoming 
very large. This indicates that the standard models for treating the rovibra- 
tional motion of diatomic molecules are no longer valid for the higher energy 
range. 

The phase shifts are the most fundamental quantity which we have com- 
puted for the scattering system (1). They shows the expected behaviors in 
terms of the variation of the angular momentum or the total collision energy. 
Our calculated total cross sections for the four individual states are similar 
to the ones obtained by one of the authors jH], and also in the similar range 
of the magnitudes to some other calculations as shown in the references jTH] 
and • However, there are still some significant differences when compared 
to all these previous study. For example, the new total cross section for 
the ground states is significantly separated from all other three dissociation 
states. In the meantime, the dissociation states all have the similar magni- 
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tudes in the total cross section. These should be more reasonable results. 

The calculated differential cross sections as demonstrated in Figures 5 
and 6 yield some useful information. For instance, for the ground state and 
with a smaller total energy, the forward scattering is a preferential one with 
some probability of backward scattering. When energy increases, however, 
the collision becomes exclusively forward and there is no chance for the head- 
on collision. This is consistent with the picture of the potential energy curve 
for the ground state. The existence of the potential energy well results in 
the quantum mechanical interference effect and therefore the splitting of the 
particle wave. With an increasing total energy or larger impact parameter, 
a classical description becomes valid for the the collision. Nevertheless, in 
contrast to the one for the ground state, the other three dissociative states 
have a dominate forward scattering. The diatomic system collides classically 
and transfers the energy to each other. This is also in accord with the shape 
of the potential curves as well as an isotropic diatomic interaction. 

Finally, we have determined the rate coefficients for the four individual 
states of the collisional hydrogen ( 2 S) and oxygen ( 3 P) system. They are all 
the increasing function of the temperature. 

In conclusion, the collisional hydrogen ( 2 S) and oxygen ( 3 P) complex is 
one of the simplest diatomic systems with many important and practical ap- 
plications. From the investigation presented in this paper, we have gained 
a better understanding of its scattering processes. Furthermore, our compu- 
tation has demonstrated how significant the shape of the potential energy 
curves is in determining the dynamic features of a molecular system. The 
range of energies being studied is also an important consideration. Therefore, 
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a dynamic theory or approach, incorporating quantum and classical mechan- 
ics transition, will be a powerful tool in the study of molecular dynamics for 
large systems. 
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Table Caption 



Table 1. The energies (in unit of eV) of the rovibrational bound states of 
diatomic molecule OH with angular momentum j = or 1 associated with 
the attractive potential well for its ground state X 2 U . 
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n,j 


E(n,j) 


n,j 


E(n,j) 


0, 


-4.349280532476708 


0, 1 


-4.344677732905585 


1,0 


-3.9007571837193575 


1, 1 


-3.8963314535284636 


2,0 


-3.4743285871919105 


2, 1 


-3.470079566811745 


3, 


-3.069862942074269 


3, 1 


-3.065790705898521 


4,0 


-2.68728756573168 


4, 1 


-2.683392721308816 


5, 


-2.326599978374491 


5, 1 


-2.3228837532003874 


6, 


-1.9878821494497058 


6, 1 


-1.9843465893508291 


7,0 


-1.6713195300265837 


7, 1 


-1.6679677223410522 


8, 


-1.3772272393753788 


8, 1 


-1.3740636477245691 


9, 


-1.1060875469879576 


9, 1 


-1.1031184700375904 


10, 


-0.8586060398134588 


10, 1 


-0.8558403676094428 


11, 


-0.6358007830952933 


11, 1 


-0.6332512865543456 


12, 


-0.43915758000136823 


12, 1 


-0.4368432791101333 


13, 


-0.27104694825436975 


13, 1 


-0.2690016425458286 


14, 


-0.1378520028115584 


14, 1 


-0.1361891922581645 


15, 


-0.05391304486808863 


15, 1 


-0.05271912534137973 


16, 


-7.919976994138954E-3 


16, 1 


-7.246986879999613E-3 
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Figure Captions 



Figure 1. The fitted potential energy curves for the four states (X 2 IT, 2 E , 4 n, 4 
of the diatomic molecule OH. 

Figure 2. The phase shifts of the s and p waves of the ground state (X 2 n) 
for the diatomic molecule OH and their variation in terms of the collision 
energy. 

Figure 3. The phase shifts of the s and p waves of the excited state 4 £~ 
for the diatomic molecule OH and their variation in terms of the collision 
energy. 

Figure 4. The total cross sections (in unit of A 2 ) for the four states 
(J\" 2 I1, 2 X - , 4 n, 4 £ _ ) of the diatomic molecule OH in the energy range from 
0.27 eV to 15 eV. 

Figure 5. The differential cross sections (in atomic units) for the ground 
state (X 2 U) of the diatomic molecule OH at the energies 0.2 eV or 3.0 eV. 

Figure 6. The differential cross sections (in atomic units) for the four states 
(X 2 II, 2 £-, 4 n, 4 £)- of the diatomic molecule OH at the energy 0.2 eV. 

Figure 7. The temperature-dependence of the reaction rate coefficient (in 
unit of cm 3 s _1 ) for the state 4 I1 of the diatomic molecule OH. 
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